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@ Process for measuring flow and determining the parameters of multilayer hydrocartwn-producing formations. 

@ The invention concerns a process tor measuring ttie flow 
as a function of time of the several layers of a subterranean 
multilayer hydrocarbon-producing formation through wtiich a 
well is drilled, and for determining the parameters of formation 
by comparison of trends in well behavior shown by actual 
measurements with trends in behavior established theoretically. 
The flow rate contributions of individual layers (1 to 5) are 
determined from cumulative measurements made above suc- 
cessive layers by a flowmeter (13) moved vertically within a 
wellbore. Experimental curves are derived that represent the 
relative variations in time of the flowrates of layers ( qj) or 
groups of layers in the formation using flowrate measurement 
points obtained when a total flowrate variation ( Q) has been 
imposed on the well at a given lime between two given constant 
values. These experimental curves then compared to theoreti- 
cal model curves (G. H) established for various values of the 
characteristic parameters of a subterranean formation, to 
determine values of the parameters of the formation. 
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Description 

PROCESS FOR MEASURING FLOW AND DETERMINING THE PARAMETERS OF MULTILAYER 
HHUOtaa HYDROCARBON-PRODUCING FORMATIONS 

r^rrsH^f^Zlsubterrlnean formations through which the wells are drilled, have long beer^ known. Although such 
productive subterranean Tormai^ns 9 considerable number of parameters characterizing 

r^rT^pi va^^o U9eS presented at the 58th Annual conference and Exposition at San Francisco. 
S?ober 5-1 1983 Ms pro^sf consists of testing each layer Individually and recording a senes of pressure 
?n^pT Si^h a orS involves at least three inconveniences. First, it takes a long time Second, the 
SetSn ^f^eTurvrs t ?^cky . there is any transfer flow between fo^^e^ion layers. Final^,. dunng 
tW,t no the well is never in an activity mode similar to a real production situation. 

AnSher method o"i^^^ -""""aver systems is to use variations of flow and pressure as a f "nctionof 

demhTa SmzeS wSr i.e . a well in which production is at a constant surface pressure and Aowra ejh s 
depth in a staomzeo we» i e . J ^ pressure at each layer for a given surface flowiate 

r pteTsTe^Te^^IlrobtredT^^^^^^ for variL successive surface fiowrates in the fom, of a 

series of pressure/flow cuives for each layer. Here, there are two incon veniences. First, not all "^f^JJ^^ 
SS«2?d^low srtuation. In addition, it was shown (Lefkovits. H.C., Hzebroek 1^.. Allen, f -E-^f Matthe««^ 
SI Ts tuTy of the Behavior of Bounded Resen/oirs Composed of Stratrfied Layers . J. Pe . Tech March 
that the^esp^ctte flowrates of the layers vaiy with time. Thus, this process is applicable only to wells 

"'fase^oS'^hr^at: o^Z S'Ls recalled, the purpose of the invention is an original process for 
deSing the chScteristic parameters of a multilayer subterranean formation. The P^^^ess consists 
21 Jn^^l Iv of detenn^ning the relative variations in time of the flowrates of layers or groups of layers of the 
S^Sto?basfd on^ow measurement points obtained when a total ftowrate variation AQ is imposed on a weH 
l^T«rS^t1 be^eenTgiven first constant value and a second given constant value, then comparing these 
SowSl vlriS to me behavior of a theoretical model established for various values of the charartenstic 
SmeteTs S a suSe^anean formation and deducing the values of the parameters of t^e formation invo^^^ 
from^hose associated with the behavior of the theoretical model which best coincide with the expenmental 

^'°Such J D°SJess makes it possible to determine the characteristic parameters of a subterranean multilay r 
formatton S S the known fact that variations of the respective flowrates of the layers of the 
forSon a?e, dur^^^^^^^ period immediately following the well flowrate change, sensitive to wall skin and 

Tth^rcra^s^a^^^^^^^^^^^^ 

flowmeter has been lowered. 

FiQure 2 represents a curve obtained by moving the flowmeter in the weH. , . „ 

PigSes 3 and 4 si^L two sets of flow cun/es prepared based on curves such as those m Figure 2. 
F gure 5 represents an experimental well pressure/time cun,e and the derivat^re curve. 
Figures 6 through 8 represent flowrate relative variation cun/es for layers relative to totaJwell flow^ 
distinct groups of ?ayers or zones relative to total well flow and layers relative to the total flowrate of the 

FigT/l'sh^ifan on S?o5SinTfSrm2- containing save., oil-bearing layers i.e.. five layers 1. 2^ 
A fnrt ^ ThP intermediate lavers 12 34 45 separating layers 1 and 2. 3 and 4. and 4 and 5 respectively, have a 
'ce^irn ve^a preS?^^^^^^^^^ hat there may be o'il flow through these intemiediate layer.. On the other 
hand lavlr 23 beLen layers 2 and 3 is impermeable and there is not oil flow between these two layers. Each 
5roi oftayers be^een which vertical oil flow can occur and which is isolated by impemrieable layers is called 
f "zone" In this example, the formation includes two zones Z1 and 22, Z1 being composed of layers 1 and 2 
andS of lave s I Sd 5. By definition, there cannot be any vertical transfer of oil behjveen two zones Tj^^ 
d?vtsSi oMhe subterranean fomiation into layers and zones has proven quite advantageous for the 
Spretation of the results obtained and is one of the characteristics of this invention. The various ayers ^d 
Snes can be ?denWied by making a recording of the flowrate as a function of depth through the formation. 

"wTn rw1» !:%:1:'::;:L^o:T.^^^^^^^^ on f,ow to the surtac through Its production string 
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11. The annular space between the casing 10 and th production string 11 is sealed off by the packer 9. The 
partial flowrat s q1 to q5 of layers 1 through 5 make up th total flowrat . Th total ftowrate Q measur d abov 
the formation is the sum of flowrat s q1 to qS. It should bo noted that the fl wrate measured on the surface 
may differ due to the wellbore storage effect. These flowrates are identical if the effect is zero. 

The process according to the invention uses the variations in time of the pr ssure p in the well and of the 5 
partial flowrates q1 to qS of the various layers resulting from a modification made to the total well ftowrate Q. 

The flow measurements are made using a flowmeter 13 (for example, as described in French patent No. 74 / 
22 391) lowered into the well at the end of a cable 14. then moved vertically several times in a sweeping 
movement throughout the total depth of the formation. When it is immediately above a layer, the flowmeter 
measures the cumulative f lowrate of that layer and those below it. With each passage, a curve such as the one W 
shown in Figure 2 is recorded, indicating the flowrate measured as a function of depth, from which the 
cumulative flowrates of the various layers 1 through 5. i.e., q5. q5 + q4. qS + q4 + q3. etc., recorded in front 
of intermediate layers 45. 34, 23. etc.. can be deduced. The times at which these flow measurements are made 
are also recorded. This makes it possible to trace the curves representing the variations in cumulative flow as a 
function of time. Figure 3 is a semi-logarithmic representation of such a set of curves, with the flowrates 75 
expressed in barrels per day (1 barret 158.98 litres). Using these curves and simple subtraction, it is 
possible to trace the set of curves in Figure 4. which represent the variations of the flowrate specific to each 
layer 1 through 5 as a function of time. 

In Figures 3 and 4. the measurement times seem to be the same for all the layers. In fact, because of the 
sweeping movement of the flowmeter 13, these points are offset from one cun/e to the next. This obviously has 20 
no effect on the curve tracing operations. 

The cable used may or may not be electrical. If it is, the data from the flowmeter are transmitted through the 
cable to the surface to be recorded and processed. When the cable used is a simple "piano string", the data 
are recorded by a downhole recorder with memories. Such a recorder is, as an example, described in British 
patent application No. 82 31560. -255 

It may be helpful to use several flowmeters connected end-to-end so as to record the respective flowrates 
of several layers at once or that of a single layer at very short intervals of time. 

Pressure measurements are made using a pressure gauge (for example, as described in the French patent 
published under the No. 2 496 884) which can be installed stationary either at the well-head or (as represented 
at 16) at the top of the formation or be connected to the flowmeter 13 (at 16') . In the latter case, it must be kept 30 
in mind that the pressure gauge is subjected to the pressure of an oil column of variable height. Measurements 
of the pressure p in the well as a function of time are obtained in this way. Like the flowrate measurements, the 
pressure measurements are transmitted by electrical cable to the surface to be recorded or are recorded in 
the well using a recorder. 

In its measurement phase, the process according to the invention consists essentially of varying the well 55 
flowrate Q by a quantity AQ at a time t1 . and measuring the pressure and the flowrates of the respective layers 
of the formation just prior to time t1, then for a certain period thereafter. The measurements taken after time t1 
make it possible to prepare the sets of flow curves in Figures 3 and 4 and the pressure curve in Figure 5 as a 
function of the time At lapsed after time t1. More precisely. Figure 5 represents the variations of the quantity 
Ap x (Q/AQ) with Ap designating the pressure differential measured between times t1 and t1 + At. The 40 
pressure scale is expressed in psi (1 psi = 6.9 kPa approx.). The variations of the derivative of the aforesaid 
quantity Ap are also represented. The abscissa and ordinate scales are logarithmic. 

The table at the end of the description gives an example of simulated values of variations of pressure Ap (in 
psi) as a function of time At (At being expressed in days and counted from time t1 ) . The values of the variations 
of the derivative (Ap)' are also indicated, calculated as explained below, as are the flowrate measurements q1 45 
to q5 (expressed in barrels/day and represented on Figure 4) of the five layers of the formation. 

The derivative (Ap)' is calculated as a function of the logarithm of At. i.e.: 



d (log ^t) 

The method of calculation and the interpretation of the derivative data are described in the published French 55 
patent application No. 83 07075 dated 22 April 1983. 

The negative values obtained for Ap and (Ap)' can be explained by the overlay principle which is well known 
to specialists. Briefly, in order for the measurements to be usable and meaningful, the well flow time prior to 
the flowrate change must be very long compared to the period of time during which measurements are made 
aft r the flowrate Chang . ^ 

In order to plot the curves in Figure 5. the valves of Ap and (Ap)' from the table are multiplied by the ratio 
Q/AQ of the flowrate Q before the chang at time t1 to the variation in the flowrate AQ before and after time t1 
In the example in Figure 5. Q = 500 and AQ = 500 - 200 = 300. as the flowrate after time t1 was reduced from 
500 to 200 barrels/day. This is equivalent to normalizing the pressure values after time t1 with the values that 
would have been obtained before time t1 . The normalization of the curves is important as regards the pressure 65 
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measurements as well as the flowrate measurements because it makes possible the use of the values 
measured just before time t1 as asympototic values for very long time periods At. 

This characteristic of the invention, which is important in practice, will be explained in connection with 
Figure 7 (points PI and P2). 

For the interpretation of the experimental pressure data, a classical analysis, well known to specialists, is 
made, consisting of plotting various logarithmic and semilogarithmic graphs to diagnose the wellbore storage 
effect; the oil flow regime in the reservoir, which can be radial and considered to be infinite at the scale of the 
well; the presence of several productive layers and the presence of possible reservoir limits. Thus, in 
logarithmic scales, the wellbore storage effect is shown by a slope equal to 1 for the pressure and pressure 
derivative curves for short time periods (beginnings of curves) and the presence of a limit or boundary to the 
reservoir is shown by an increased in the pressure and pressure derivative values for long periods (ends of 
curves). These diagnostic methods are commonly used in the petroleum industry and are described, for 
example, in U.S. patent 4 328 705 and published French patent application No. 83 07075. 

The convolution of the total flow measured at the bottom of the well with the pressure can also be used to 
eliminate the wellbore storage effect from the pressure measurements (in this case, the pressure is called rate 
convolved pressure). This technique is published in "Interpretation of Pressure Build-up Test using In-situ 
Measurement of Afterflow". Journal of Petroleum Technology, January 1985. 

In its measurement interpretation phase, the process according to the invention includes the following 
parameter determination operations: 

A) kh (average product of formation permeability k x thickness h for the overall formation); 

B) kj and sj (horizontal permeability and skin coefficient of layer j, with j varying from 1 to 5 in this 
example); 

C) type and position of the external limit or boundary of the formation (which determines the extent and 
type of the formation) ; 

D) vertical permeability between layers. 

A. Determination of the parameter kh 
Based on pressure measurements, using the following formula: 



141.2 B fl 

in which: AQ is the change made in the well flowrate (expressed in barrels/day) at the time ti; 
B is the relative volume factor of the oil in the formation and at the surface (equal to the ratio between the 
volumes of oil in the formation and at the surface); 
li is the viscosity of the oil expressed in centipoises; 

(Ap)'M is the value of the derivative of the pressure p as a function of the logarithm of time in the flat part of the 
derivative curve (Figure 5). This flat portion reflects an infinite action radial flow. 

In this example, Figure 6 shows that: 
(Ap) M . (Q/AQ) = 5 
Thus, (Ap) M = 3. 
45 In addition, other measurements revealed that: 
B = 1.2 and ^ = 1. 

Thus, kR « 8.472 md.ft 
- 25.42 x 102 |im2.m. 

50 B. Determination of kj and sj 

For each layer j, the curve representing as a function of time the fraction of the variation of total flow 
attributable to the layer involved, i.e., the quantity Aqj/AQ, based on the values from the table and the curves in 
Figure 4. This results in five series of points in semilogarithmic representation (Figure 6) for the five layers 1 
through 5 of the formation, respectively, as a function of the time At lapsed after the time t1. 

55 Each series of points is then compared with a theoretical model to determine which of the curves suitably 
fits the series of points involved, at least during the initial period following the time t1 of the flowrate change. It 
has been acknowledged that during this period the model used may correspond to the absence of flow 
between layers in the formation, with an infinite external boundary. After the initial period, deviations may 
appear between the measurement points and the theoretical curve due to interlayer flow and extemal 

60 boundary effects, or overlay effects. 

The theoretical model was established based on the following formula: 
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in which ^jd is the Laplace transform of the dimensionless flowrate of layer j; 



k. = (kh)j/kh with kh 



sj is the skin coefficient of layer j; 

Ko and Ki are the modified Bessel functions of the first and second types; 
PwD is the Laplace transform of the dimensionless pressure in the well; 



while ffj 1s given by: 



In which: 



n 



J 



with: 

in which (9h)j designates the product porosity x height of layer j, n the number of layers in the formation and z 
the Laplace space variable. 

Equation ( 1 ) does not give flowrate as a function of time. To obtain it, the inverse Laplace transform given by 
the Stehfest algorithm is applied (see 'Numerical inversion of Laplace transforms", D-5. Communications of 
the ACM. January 1970, No. 1, pages 47 to 49). 

When fitting is achieved with a given curve of the theoretical model, the skin coefficient sj of the layer j 
involved, appearing in formula (1), can be deduced from it. as can Its permeability, which also appears in 
formula (1) as the product (kh)j of the permeability and height of said layer, the latter parameter being known 
by previously made log measurements, while the product ER is determined using the pressure 
measurements explained above. 

To illustrate cases which could be encountered in practice, two theoretical curves G and H (dotted tines) are 
shown in Figure 6 which do not correctly fit the series of measurement points for formation layers 1 and 2 on 
the left side of the figure, while there is a good fit on the right side (significant deviations do, however, occur at 
the extreme right due to boundary and interlayer flow effects). The examination of the position of curve G 
shows, for example, that the skin coefficient selected for it is too low and should be increased. For curve H, the 
opposite is true: the skin effect must be reduced, even though a modification of the value of the skin effect of 
curve G has an influence on the other curves. 

These operations make it possible to determine the horizontal permeability k] and the skin coefficient sj of 
each layer of the formation. 

Using a different operational mode of this invention, these parameters can be determined as explained 
below: 

It is known to specialists that the result of the mathematical operation of convolution of the derivative of the 
variations of flow q with dimensionless pressure Pd (the pression that would be obtained if no other 
parameters intervened in the formation and the well to influence the pressure value and if the flow rate were 
constant) represents th variations of the pr ssure Psi effectively measured in the well in front of the 
formation. This is expressed by the following equation: 
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T 

Psf (T) = I q" (T < t) . I>D (t) Ax 



Psf (T) being the value of the pressure variation measured in the wet! at time T. 
To obtain Pd, which is the pressure value being sought, requires the mathematical deconvoiution between 

ro the effectively measured pressure Pstand the flow. However, the results obtained by deconvolution can be 
sprinkled with significant errors if the experimental data include some noise. Convolution is thus the preferred 
operation. This is why, within the framework of this invention, the flow variations for each layer and the 
pressure variations in the well were measured. It was then shown that the convolution of the flow variations for 
each layer with the pressure variations tn the well provides the pressure response of the layer as if it were the 

75 only one producing a fluid, provided, however, that there is no interlayer flow. Thus, armed with the pressure 
response of each layer, it is possible to use the classical methods of interpretation for each, spec'ificafty the 
pressure/time curves plotted on semitogarithmic scales which make it possible to determine the permeability 
and skin effect. 

20 C, Determination of the external boundary of each zone 

The point is to determine the boundary type of each zone: 

- seemingly infinite boundary; 

- no-flow boundary, behaving like an impermeable seal, with ail the liquid flowing Into the well coming from the 
formation zone located inside this boundary; 

25 - Constant pressure boundary. 

In the first instance, it is as though there were no boundary. In the other two cases, the radius of the 
boundary must also be specified. 

For this determination, a graph [Figure 7) simitar to Figure 6 is prepared, in which each series of points 
corresponds to a zone i of the formation based on the above definition, and no longer to a given layer (of 
30 course, a zone may contain only one layer). 

In this example, there are two zones Z1 and 22 (Figure 1) and the two series of points represent, 
respectively, the following quantities: 



35 



50 



60 



65 



Aql + Aq2 Aq3 + Aq4 + ZqS 

« 



as a function of the time At. The values from the table are used to plot the experimental curves In Figure 7. 
40 In addition, theoretical models corresponding to the above formula (1). as well as the following formulas, 
were used: 

-fS «IjD - — — (2) 



and 



55 %-^^j''j*iOj 



In these formulas, formula (2) being already known, the quantities are defined as follows: 



BNSDOCID: <EP 0217664AlJ„> 



* 0217684 

♦01K<2> = ^C«'k>Il(^eD''k> * Io<''k>'^l^reD'^k> 



70 



in which lo, li, Ko and Ki are modified Bessei functions of the first and second type and reo >s the 
dimensionless external radius of the formation. 

Formula (1) refers to the case of a boundary behaving as H it were infinite, formula (2) a no-flow boundary 
and formula (3) a constant pressure boundary. '5 

Figure 7 shows the fitting achieved (in the initial period) between the series of points corresponding to 
zones 21 and 22 and the curves defined based on formula (2). It can be concluded that the boundary of the 
zones under study is of the "no-flow" type. 

If the fitting is achieved with a curve ending with a horizontal, such as curve K con-esponding to formula (3). 
sketched at the top right of Figure 7. the boundary is of the "constant pressure" type. If the curve ends with a 20 
slight downward bend, like cun/e L resulting from formula (1), the boundary is of the apparently infinite type. 

Indeed, in these operations relative to the various zones of the fomiation, there is no need to envision a 
vertical transfer flow situation, since such transfers are by definition inexistent between zones. 

In the case represented in Figure 7. showing a no-flow boundary, i.e., a situation in which the production 
volume of the formation is limited, the straight portion of the cun/e for each zone tends towards a value equal 25 
to the product <ph for the zone, i.e., in this example, 0.4 for zone Z1 and 0.6 for zone Z2. These values are in 
addition known due to previous logging operations. 

Figure 7 also shows that the straight portion of the curves deviates from the experiinental points. This is 
because of an overlay effect due to the fact that the well in question was placed in production for 200 hours 
(8.33 days) and then its production flow rate was reduced (from 500 to 200 barrels/day) at time t1 , for another 30 
200 hours. However, if a measurement is made at the end of the first 200-hours period, just prior to time t1 . the 
results obtained can be recorded on the figure (points P1 and 92) and considered as measurement points 
obtained after the initial period, after time t1 . without overiay effect. As can be seen, the theoretical curves are 
very close to these points. 

It results from the foregoing remark that, in practice, there is no need to make measurements at times 35 
distant from time t1 . since measurements made Just prior tc that time can advantageously replace them. Thus, 
the measurements at the points located approximately between 10° and 10^ days after t1 need not be made, 
and this considerably shortens the total time required for well measurements. 

As a result, it is possible to determine using only the measurements taken just before t1 whether the 
boundary is of the no-flow type (if the measurement points con-espond to the known values of 0h) or not.(if the 40 
measurement points do not con^espond). since this effect depends solely on boundary conditions. 

If the boundary is recognized to be non-infinite, its radius is determined by finding the radius value which 
leads to the best fit between the model curves and the measurement curves in the period following the Initial 
period. The measurement points recorded just prior to the change In the weirs flow rate are also very useful in 
this phase of determination of formation parameters. 

As in the case of the determination of permeability and skin effect, there is a second possible method using 
convolution operations. It has, in fact, been shown that the convolution of the flow variations of each zone with 
the pressure variations in the well provides the pressure response of the zone involved. As in the case of the 
individual layers, this harks back to a classical well test interpretation, particularly for the determination of the 
boundary of each zone. ^ 

D. Determination of interlayer permeability 

Having determined the horizontal permeability and the skin coefficient of each layer of the formation, the 

type and location of the zone boundaries, the vertical permeability between layers remains to be determined. 

This is done by means of an analysis in each zone of the formation of the flow rates of the layers of the zone as 55 

compared to the zone's total flow rate- 
As is shown in Figure 8. we show as a function of At. still using a semilogarithmic representation and based 

on the table data, the quantities Aq q!/AQi in which AO* designates the flow rate variation of zone i and Aqij the 

flow rate variation of layer j belonging to zone i. Figure 8. as an example, is limited to the measurements for 

zone Z1. composed of layers 1 and 2. the values for which are indicated in the table (page 19). 60 
The th oreticat model used, established for the case in which there are transfer flows between layers, is 

derived from the following formula: 

65 
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m. 

1 



in which: 

10 QjD is the dimensionless flow rate of layer j of zone i. which contains mi layers; 
Co is the dimensionless wellbore storage constant; 

Kj - ratio of the product permeability x height for layer j to the average product Eli permability x height for 

zone i: 

ajci represents the quantity ak for zone i; 
15 Ok are the roots of the equation yn ^ 0 in which yj »s a polynomial defined by recurrence by: 

20 for j = 2 n with yo = 1 and yi = an. ajk designating the elements of the matrix ajkl : 



for k==j-l ; j>l, 

2 

k.<T -«.z-X. T-X-f for k=j, 

X^, ^ ^"-^ ^ for k=.j+l ; j<n 

I for k9£j-l,j, or j+1 



a is a coefficient relative to layer j, root k, for zone i. defined by the formula: 



30 



k.=l 

35 1 

bki Is an external boundary condition coefficient defined by the formulas: 
b»ci - 0 for an apparently infinite boundary 



= KiC«^ki^eD>^^l<:^ki':eD> 

for a no-flow boundary 



for a constant pressure boundary, 

while B I* is related to A j[* by the equation: 

50 = b 

A 1^ being determined based on well conditions. 

The skin effect coefficient values obtained in interpretation phase B) are used, keeping In mind the type and 
55 location of the formation's external boundary as determined in phase C). Finally, a set of values is sought for 
the parameters X\ of interlayer permeability between layers j and j + 1 of each zone i such as to achieve good 
fit of the curves for all the Aqj / AQiratios considered. 

More precisely, it can be noted that the appearance of the curves in the left-hand portion of the figure 
depens on permeability and skin effect, while the right-hand side of Figure 8 depens also on the type of 
60 boundary and transfer flows. Since the permeability, skin effect and boundary type are known, the only 
remaining parameter is transfer flow, for which different values are tried until a good curve fit is achieved. 

These operations are repeated for each of the zones of the formation, so as to determine the transfer flow 
parameters for all layers. 

The set of calculations and curve fitting op rations just described as part of the process according to the 
6b invention can be done my hand or, pref rably, by a digital calculator. In the first instance, sets of typical curves 
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are traced using the equations given above. These sets of curves are a graphic representation of the behavior 
of the theor tical models. A digital catculator can aiso be used to select the values of the paranneters being 
sought which corr spond to a perfect fit between the the r tical and experimental variati ns of the various 
functions of the pressure and flow rates (variation of pressure, of the derivative of pressure, of the fraction of 
the variation of the total flow rate for a given layer and for a giv n zor>o and of the fraction of th variation of flow 
rate of a layer as compared to the flow rate of the zone to which It belongs, all as a function of time) . 
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1. Process for the measurement of the specific flowrate in a well, as a function of time, of the various 
layers of a subterranean multi-layer hydrocarbon-producing formation, characterized by the fact that it 50 
consists of lowering a flowmeter (13) into the well (10), moving it successively before each of the layers (1 
through 5) of the formation and recording the flowrate values It provides when located immediately above 
each layer and the measurement times, and then deducing the flowrate specific to each layer (1 through 

5) from the cumulative flowrates thus recorded. 

2. Process according to claim 1. characterized by the fact that the flowmeter (13) is repeatedly raised 55 
and lowered throughout the entire depth of the formation, so as to determine the flow variations of each 
layer as a function of time. 

3. Process according to claim 1 or 2, characterized by the fact that a pressure gauge (16) Is connected 
to the flowmeter (13) to provide pressure measurements in the well (10) which are concomitant with the 

flow measurennents. 60 

4. Process for determining the characteristic paramet rs of a multilay r subterranean hydrocarbon- 
producing formation thr ugh which a well is drilled, 

characterized by the fact that it consists essentially of det rmining the relative variations in time of the 
flow rates of layers or groups of layers of the formation with r spect to flow rate measurement points 
obtained when a total flow rate variation AQ has been imposed on the well at a time t1 between a first 65 
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given constant value and the second given constant value, then comparing those flow rate variations with 
the behavior of a theor tical model establish d for various values of the characteristic parameters of a 
subterranean formation and deducing the values of the parameters of the formation involved from those 
associated with the behavior of the theoretical model which best fit the experimental flow rate variations. 

5. Process according to claim 4, characterized by the fact that well pressure values measured at the 
same time as the flow rates of the various layers are also used, making it possible to calculate the product 

permeability x formation thickness. 

6. Process according to claim 4 or 5 characterized by the fact that for each layer j of the formation the 
variations of the fraction Aqj /AQ representing the ratio of flow rate variations Aqj of layer j to the variation 
AQ of the total well flow rate Q are determined as a function of time At, and the horizontal pemneability 
and skin coefficient parameters kj and sj of the layer are deduced by comparing said variations Aqj/AQ 
with the behavior of a theoretical model. 

7. Process according to claim 5 characterized by the fact that the variations of flow rate of each layer j 
are convoluted with the well pressure variations, which makes it possible to obtain the pressure response 
of the given layer as if it were alone, from which the permeability and the skin effect coefficient of the layer 
can be deduced using classical methods of Interpretation. 

8. Process according to any of claims 4 through 7 characterized by the fact that for each zone i of the 
formation, i.e., for each group of productive layers contained between two impermeable Intermediate 
layers, a determination is made of the variations as a function of time, of the fraction lAcf I AQ 
representing the ratio of the flow rate variations of said zone 1 to the variation in the total flow rate Q of the 
well, and the type and position of the external boundary of said zone are deduced from the comparison 
between said variations of said fraction and the behavior of a theoretical model. 

9. Process according to any of the claims 4 through 7 characterized by the fact that the flow rate 
variations of each zone i are convolved with the well pressure variations to obtain the pressure response 
of the zone involved, from which the type and position of the extemal boundary of said zone can be 
deduced by classical methods of interpretation. 

10. Process according to any of the claims 4 through 9 characterized by the fact that for each layer j of a 
zone i of the formation, i.e.. of a group of productive layers contained between two impermeable 
intermediate layers, a determination is made of the variations of the fraction Aqj/AQTepresenting the ratk) 
of the flow rate variations of layer j to the flow rate variation of zone i, and the interiayer permeability 
parameters of the zone involved are deduced by comparison between the said variations of the said 
fraction and the behavior of a theoretical model, for each of the zones in the formation. 

11. Process according to any of claims 6 through 10 characterized by the fact that the curves 
representing the variations of said fractions are plotted and then fitted to the curves representing said 
theoretical model of the characteristics of the subterranean formation. 

1 2. Process according to any of claims 4 th rough 1 1 characterized by the fact that instead of the flow rate 
measurements made after a relatively long lapse of time A t after the time t1 when the well flow rate is 
changed, the measurements made just prior to time t1 are used, where the well has been flowing for a 
period essentially equivalent to the aforesaid time lapse. 

13. Process according to any of claims 4 through 12 characterized by the fact that the measurement of 
the specific flow rate as a function of time of the various layers of the formation is achieved by lowering a 
flowmeter (13) into the well (10) and moving it successively in front of each of the layers (1 through 5) of 
the formation and recording the flow rate values and time readings it provides when it is located 
immediately above each layer, and then deducing the flow rate specific to each layer (1 through 5) from 
the cumulative flow rates thus recorded. 

14. Process according to claim 10 characterized by the fact that the flowmeter (13) is repeatedly reused 
and lowered throughout the depth of the formation. 

15. Process according to claim 5 and claim 13 or 14 characterized by the fact that a pressure gauge (16) 
is added to the flowmeter (13) to provide well presure measurements concomitant with the flow rate 
measurements. 

16. Process according to claim 15 characterized by the fact that said pressure gauge is in afixed position 
in the well. 



10 



•6401653 0 




J. a 



CumulaNVe 
Q Ftowrate 




BNSDOCID: <EP 02176S4A1_L> 




BNSDCX^ID: <EP 0217684A1 J_> 



3/4 



86401653 0 




BNSDOCID: <EP 02176S4A1J_> 




BNSDOCID: <EP 0217684A1 J_> 



European Patent 
Office 



EUROPEAN SEARCH REPORT 



Application number 



EP 86 40 1653 



DOCUMENTS CONSIDERED TO BE RELEVANT 




Cmgofy 


Citation ol doeumant wtttilndiCBtion. wtwn appropriaH. 
of ralavarvt paasaQaa 


Relevant 

todatm 


CLASSIFICATION OF THE 
APPLICATION (Int CI.' ) 


X 


FR-A-2 434 923 (SOCIETE DE 
PROSPECT I ON ELECTRIQUE 
SCHLUMBERGER) 

* Page 6, lines 12-30; page 8, 
ixnes / — lo; page ixnes -Zo; 
page 2, lines 4-14; page 3, lines 
14-21 * 


1-15 


E 21 B 47/10 
E 21 B 47/06 
E 21 B 49/00 
E 21 B 49/08 


Y 




16 






FR-A-2 544 790 (ETUDES ET 

FABRI CAT I ONS F LOPETROLj ) 

* Page 9, line 26 - page 10, line 

5 * 


16 




A 




4-15 










TECHNICAL FIELDS 
SEARCHED (Int. a «| 


X 


FR-A-2 238 836 (SOCIETE DE 

PROSPECT I ON ELECTRIQUE 

SCHLUMBERGER) 

* Page 4, lines 5-19 * 


1,2 


E 21 B 


Y 




3 




A 




4-16 




Y 


US-A-4 015 194 (EPLING) 

* Column 4, lines 13-21; column 

3, lines 5-11 * 


3 




Th» pratwit March raport liM bMO drawn up (or •« daliRS 







THE I^GUE 



Oaiaof . _ . 
04-11 



Examiner 

SOGNO M.G. 



CATEGORY OF OTEO DOCUMENTS 

X : particularly relevant if talcan alona 

Y : particularty relevant if comtMned with another 

document f the same category 
A : technological backgrourHl 
O : non-written disclosure 
P : imermediate document 



T : theory or principle urMlerlyif^ the invention 
E : eartlsr patent document. tMJt published n. or 

after the filing date 
D: document cited in the application 
L : document cited for other reasons 

& : member f the same patent family, corresponding 
- document 



BNSDOCID: <EP 0217684A1 J_> 




European Patent 
Office 



EUROPEAN SEARCH REPORT 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Application number 

EP 86 40 1653 
Page 2 



Category 



Citation of document with indication, whefe appropriate, 
of relevant F 



Relevant 

to claim 



CLASSIFlCATIOfI OF THE 
APPLICATION <lnt. a.4) 



P,X 



FR-A-2 101 037 (SGCIETE DE 
PROSPECT I ON ELECTRIQUE 
SCHLUMBERGER) 

* Page 6, line 18 - page 7, line 
11 * 



EP-A-0 176 410 (SGCIETE DE 

PROSPECTION ELECTRIQUE 

SCHLUMBERGER) 

* Whole document * 



1.2 



1-16 



7, 



A,D 



PETROLEUM ENGINEER 
INTERNATIONAL, vol. 51, no. 
June 1979, pages 129,130, 
Dallas, Texas, US; E.H. 
TIMMERMAN: "Various methods used 
to average permeability" 
* Whole document * 



US-A-4 328 705 (GRINGARTEN) 
* Whole document * 



1,4-16 



The present search report has been drawn up for all claims 



TECHWCAL FIELDS 
SEARCHED <lnt. O*) 



4-16 



Place of search 

THE HAGUE 



Date of completion of the eearch 

04-11-1986 



Examiner 
SOGNO M.G. 



CATEGORY OF CITED DOCUMENTS 

X : particularly relevant if taken atom 

Y : particularly relevant if comt>ir>ed with another 

document of the same cateooty 
A : technological background 
O : non-written disclosure 
P : int rn>ediate document 



T : theory r principle underlying the invention 
E : eartier patent document, but published on. or 

after the filing dat 
D : docuntent cited in the application 
L : document cited for otfier reasons 

& : member of the same patent family, corresponding 
document 



BNSDOCID: <EP 02176e4A1J_> 



